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' ABSTRACT 

This  report  describes  the  design  of  a compact  plenum  chamber 
the  function  of  which  is  to  receive  air  from  a rectangular  duct  and 
deliver  it  to  an  annular  space  with  low  energy  losses.  Turning  vanes 
formed  from  a single  thickness  of  sheet  steel  ensure  that  air  leaves 
the  chamber  essentially  perpendicular  to  the  plane  of  the  annulus.  A 
FORTRAN  program  is  presented  which  plots  the  two  dimensional  develop- 
i ment  of  the  turning  vane. 
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1.  INTRODUCTION 

This  report  describes  the  design  of  an  air  supply  manifold 
for  the  film  cooling  facility  at  the  Defence  Research  Establishment 
Suf field  (DRES).  The  manifold  was  required  to  receive  air  from  a 
rectangular  duct  and  deliver  it  to  an  annular  space  with  lower  energy 
losses  than  had  been  experienced  with  previous  manifolds.  Ideally, 
velocity  vectors  would  be  perpendicular  to  the  plane  of  the  annulus 
and  uniform  in  magnitude;  however,  devices  downstream  of  the  annulus 
could  be  used  to  smooth  out  the  velocity  profile.  The  essence  of  the 
problem  was  to  move  air  efficiently  through  what  added  up  to  be  a 90 
degree  turn. 

Figure  1 shows  the  DRES  Film  Cooling  Facility.  The  Orenda  8 
turbojet  engine  draws  air  through  an  intake  on  the  south  side  of  the 
building  and  expels  hot  exhaust  gases  into  the  exhaust  duct  on  the  north 
side  of  the  building.  Excess  gas  is  diverted  by  the  continuously  adjust- 
able bypass  unit.  After  the  bypass  the  exhaust  flows  through  24  feet  of 
2-foot  inside  diameter  pipe  to  the  film  cooled  test  section  which  is 
approximately  2 feet  in  diameter  by  2 feet  in  length.  The  15  HP  centri- 
fugal fan  drives  air  through  approximately  40  feet  of  20-inch  diameter 
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duct,  and  the  manifold,  to  the  test  section.  Convective  heat  transfer 
from  the  hot  exhaust  gas  to  the  walls  of  the  test  section  is  reduced  by 
the  film  of  cool  air. 

Figure  2 shows  the  original  manifold,  called  Manifold  I, 
with  tangential  and  radial  air  entry,  and  the  improved  Manifold  II. 

The  film  cooled  test  section  and  cooling  air  duct  are  shown  in  phantom 
lines.  It  was  found  that  systems  using  Manifold  I had  unsatisfactory 
total  cooling  air  flow  and  velocity  distribution  in  the  test  section. 
Efficiency  of  the  air  supply  system  was  defined  as  where 

is  velocity  pressure  at  the  annulus  and  P is  the  fan  total  pressure. 

rT 

This  ratio  was  experimentally  determined  to  be  4 - 28%  for  various 
systems  using  Manifold  I.  Data  revealed  that  Manifold  I was  the  largest 
source  of  inefficiency  in  the  air  supply  system.  Redesign  of  the  mani- 
fold had  the  greatest  potential  for  improvement  of  efficiency. 

A lower  entry  velocity  and  more  streamlined  interior  would 
reduce  turbulent  losses.  Tangential  velocity  components  at  the  annular 
exit  of  Manifold  I were  not  good  for  film  cooling  and  it  would  be 
advantageous  to  convert  them  into  axial  velocity  components.  These 
guidelines  were  followed  in  the  design  of  Manifold  II. 

2.  APPROACH 

Two  approaches  were  considered.  First,  a large  plenum  chamber 
could  be  used.  The  path  of  the  air  from  entry  to  outlet  of  the  chamber 
would  not  be  known  precisely  but,  because  of  low  velocities,  low  resist- 
ance and  good  distribution  could  be  expected.  Second,  the  manifold 
could  be  compact  and  direct  relatively  high  speed  flows  over  stream- 
lined surfaces.  It  was  not  known  what  size  of  large  plenum  chamber  was 
necessary.  Support  posts  for  the  exhaust  duct  in  the  film  cooling 
facility  would  interfere  with  any  large  chamber.  Therefore,  it  was  de- 
cided to  build  a compact,  streamlined  manifold. 
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3.  FEATURES 

Figure  3 shows  the  general  arrangement  of  Manifold  II.  Further 
detail  is  contained  in  Figures  4 and  5.  The  inner  shell  is  made  from 
sheet  steel  to  withstand  high  exhaust  gas  temperatures.  The  outer  shell 
is  made  from  plywood  for  ease  of  fabrication.  Temperatures  occurring  in 
the  outer  shell  are  quite  low.  A layer  of  fiberglas  over  the  plywood 
adds  strength. 

The  cross-sectional  area  of  the  air  inlet  is  144  square  inches, 
compared  to  64  square  inches  for  Manifold  I tangential  entry  and  91 
square  inches  for  Manifold  I radial  entry.  The  air  inlet  velocity  to 
Manifold  II  is  therefore  lower  than  to  either  of  tho  Manifold  I variations. 

Between  inner  and  outer  shells  of  Manifold  II  is  a trapezoidal 
passage  (Figure  5)  the  cross-sectional  area  of  which  decreases  proportion- 
ally to  the  angle  travelled  around  the  manifold.  Constant  velocity  is 
maintained  in  the  trapezoidal  passage  by  decreasing  area  as  air  bleeds  off 
to  the  turning  vanes  and  test  section.  The  function  of  the  turning  vanes 
is  to  eliminate  tangential  velocity  components.  All  the  vanes  are 
identical  in  design  because  velocity  in  the  trapezoidal  passage  is  con- 
stant. 

The  inner  shell  shown  in  Figure  5 is  conical  so  that  most  of 
the  manifold  is  separated  from  the  exhaust  pipe  by  fiberglas  insulation 
and/or  air  space.  Width  of  the  turning  vane  trailing  edge  is  fixed  by 
choice  of  the  maximum  slot  width  for  the  film  cooled  test  section.  Pre- 
liminary calculations  showed  that  if  the  width  of  the  vane  was  the  same 
at  the  leading  edge  as  at  the  trailing  edge,  the  effective  cross-section- 
al area  between  the  vanes  at  the  trailing  edge  would  be  3.28  times  that 
at  the  leading  edge.  When  the  walls  of  a diffuser  diverge  too  quickly, 
air  flow  tends  to  separate  from  the  walls.  Flow  separation  will  also 
occur  on  the  inner  surface  of  an  elbow.  The  passage  between  vanes,  as 
shown  in  Figure  3,  is  both  a diffuser  and  and  elbow.  Thus  the  ideal  uni- 
form velocity  profile  at  the  exit  from  the  vane  passage  is  unlikely  to  be 
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achieved.  There  will  be  a region  of  low  velocity  near  the  "Inner" 
vane  and  a region  of  high  velocity  near  the  "outer"  vane.  Increasing 
the  width  of  the  vane  leading  edge  would  tend  to  relieve  the  velocity 
profile  problem  because  It  would  make  the  passage  less  divergent.  How- 
ever, another  factor  complicates  the  situation.  Figure  5 shows  that 
after  leaving  the  vane  passage,  the  air  must  make  another  slight  turn  to 
flow  parallel  to  the  test  section  wall.  Increasing  the  width  of  the 
vane  leading  edge  would  increase  the  average  angle  through  which  the  air 
would  have  to  turn  when  leaving  the  vane  passage.  A compromise  leading 
edge  width  was  chosen  so  that  the  previously  mentioned  area  ratio  be- 
came 2.68  and  the  turning  angle  at  vane  exit  was  19  degrees. 


4.  TURNING  VANES 

Airfoil  shaped  vanes  would  be  most  desirable  but  casting  or 
forging  facilities  were  not  readily  available.  Vanes  were  bent  from 
single  thicknesses  of  sheet  steel. 

Cylindrical  and  rectilinear  coordinates  are  defined  in  Figure  6. 
The  theoretical  surface  of  the  vane  is  generated  by  a radial  line  seg- 
ment, the  generatrix,  simultaneously  rotating  in  the  9 direction  and 
translating  in  the  z-direction.  Intersections  with  conical  surfaces  form 
the  inner  and  outer  edges  of  the  vane. 


Energy  losses  due  to  turbulence  are  minimized  by  ensuring  that 
the  tangent  plane  to  the  vane  leading  edge  is  coplanar  with  the  velocity 
vector  of  the  air  leaving  the  trapezoidal  passage  (Figure  7).  V is  the 
velocity  vector  of  the  air  at  the  vane  leading  edge.  V , V and  V are 
its  tangential,  axial,  and  radial  components  respectively.  In  order  that 
the  leading  edge  tangent  plane  be  coplanar  with  V,  the  angle  e must  equal 


tan"^ 


V 

a 


Figure  8 illustrates  some  of  the  quantities  used  in  the  cal- 
culation of  V and  V . Volume  flow  rate  into  the  manifold  through  Aq 
a t 

equals  volume  flow  rate  out  of  the  manifold  through  Ai  assuming  constant 
air  density  and  no  leaks. 
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7r(r  2 _ ^2) 

10  1 

By  the  principle  of  conservation  of  angular  momentum,  the  pro- 
duct of  radius  times  tangential  velocity  is  constant,  neglecting  change 
in  momentum  due  to  friction  and  turbulence  losses;  hence  at  the 
mean  radius  of  the  vane  leading  edge  is  estimated  as 


Thus 


Note  that  B depends  only  on  the  geometry  of  the  manifold  and  does  not 
vary  with  flow  rate.  Using  average  values  of  radius  and  velocity  results 
in  small  errors  for  points  on  the  vane  leading  edge  other  than  the  mid- 
point. 

The  orientation  of  the  leading  edge  tangent  plane  is  now  known, 
and  the  trailing  edge  tangent  plane  must  be  coplanar  with  the  z axis.  The 
scheme  chosen  to  describe  the  curved  surface  joining  the  two  tangent  planes 
is  a curve  of  z versus  the  r9  product  as  constructed  in  Figure  9. 

p = (>^1  <^10)  . Prom  the  curve  is  obtained  a set  of  (9,  z)  pairs. 

2 

Next,  r values  are  attached  to  the  (9,  z)  pairs  for  points  on 
the  inner  and  outer  edges  of  the  vane.  These  edges  are  the  intersection 
of  the  vane  with  conical  surfaces.  See  Figure  10.  For  the  inner  edge. 


and  for  the  outer  edge 


r 


As  explained  by  Reference  1,  the  vane  is  a warped  surface  which 
cannot  be  developed  on  a plane.  It  is  necessary  to  approximate  the  warped 


1 


1 


J 
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surface  with  a series  of  triangles  so  the  vane  can  be  laid  out  on  a 
flat  sheet.  Figure  11  illustrates  the  three-dimensional  vane  and  its 
two-dimensional  pattern.  Points  1 to  18  lie  on  the  theoretical  inner 
and  outer  edges  of  the  vane  and  their  (r,  0,  z)  coordinates  are  known. 
The  lengths  of  all  the  line  segments  are  claculated  using  the  Pytha- 
gorean theorem  for  three  dimensions. 

The  flat  layout  is  started  with  segment  1-10  along  the 
y-axis.  Arcs  having  the  lenths  of  segments  1-11  and  10-11  could  be 
struck  from  points  1 and  10  respectively  to  locate  point  11.  Alter- 
natively, the  angle  between  segments  1-10  and  1-11  could  be  calculated 
using  the  law  of  cosines,  and  then  point  11  located  trigonometrically. 
The  latter  method  lends  itself  to  computer  programming.  Point  2 is 
located  after  point  11,  then  point  12  and  so  on. 

A quadrilateral  9-19-20-18  is  added  to  the  vane  coplanar 
with  the  z-axis  to  help  establish  flow  in  that  direction.  Point  19 
is  separated  1/4  inch  from  the  inner  cone  to  reduce  heat  transfer  and 
problems  due  to  thermal  expansion. 

A FORTRAN  program  was  written  following  the  methods  outlined 
above  to  calculate  coordinates  and  plot  a two-dimensional  vane  pattern. 
It  requires  as  input  the  (9,  z)  coordinates  of  points  1 to  18  in- 
clusive, and  the  set  of  point-number  pairs  for  the  end  points  of  line 
segments.  The  variable  definitions  and  program  listing  are  included 
as  Appendix  A,  while  the  full  scale  vane  pattern  is  shown  in  Figure  12. 
Note  that  the  pattern  has  been  reflected  about  the  x-axis  because  the 
tangential  velocity  in  the  actual  manifold  is  opposite  to  that  imagined 
for  the  development  of  the  equations  describing  the  vane.  Another 
FORTRAN  program,  not  included  here,  plotted  the  lines  of  contact  of  the 
vane  on  the  inner  and  outer  conical  surfaces.  These  lines  aided  con- 
struction of  a three-dimensional  vane  model,  and  assembly  of  the  vanes 
on  the  conical  surfaces. 
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5.  CONCLUSION 


P 

Efficiency  ^ of  the  cooling  air  system  using  Manifold  II 
was  0.52  compared  to  O^B  for  systems  using  Manifold  I.  Thus  the  major 
objective  of  the  design  was  achieved.  The  velocity  profile  had  irregu- 
larities in  direction  and  magnitude,  but,  as  previously  mentioned,  de- 
vices downstream  of  the  manifold  tended  to  smooth  out  the  profile. 
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FORTRAN  PROGRAM  TO  PLOT  TURNING  VANE  PAHERN 


Variable 

Definition 

Mode 

Class 

R(N) 

radius  r to  point  N (inches) 

real 

1 

TTA(N) 

angle  0 to  point  N (degrees) 

tf 

II 

Z(N) 

height  z of  point  N (inches) 

II 

II 

ALONG  (L) 

length  of  line  segment  L (inches) 

11 

II 

X(N) 

x-coordinate  of  point  N in  plot 

II 

II 

Y(N) 

y-coordinate  of  point  N in  plot 

II 

II 

N 

point  index  number 

integer 

0 

L 

line  segment  index  number 

II 

It 

N1 

point  number  at  beginning  of  line  segment  L 

11 

II 

N2 

point  number  at  end  of  line  segment  L 

II 

II 

ALPHA 

angle  of  line  segment  L from  reference  line 
1-10  (degrees) 

real 

II 

J 

vertex  angle  index  number 

integer 

II 

C 

cosine  of  vertex  angle  V 

real 

It 

V 

vertex  angle  between  two  line  segments 
(degrees) 

II 

II 

B 

variable  to  indicate  whether  J is  even  or  odd 

II 

II 

AB 

smaller  angle  between  line  segment  9-18  and 
y-axis  (degrees) 

II 

II 

I 

index  number  used  to  set  order  of  plotting 
of  line  segments 

i nteger 

II 

D 

variable  to  indicate  whether  I is  even  or  odd 

real 

II 

XX 

x-coordinate  of  final  pen  position 

II 

II 
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TURNING  VANF  COMPOSED  OF  TRIANGULAR  ELEMENTS 
DRFSMES  D.A.  RAYLV 

THIRD  VERSION 

FULL  SCALE  PATTERN  OF  MANIFOLD  II  VANES 
DIMFNSION  R ( 18  ) .TTA(  18)  iZ  ( 18  ) tALONGOS  ) *X(  20)  »Y  <20  ) 


READING  SET  OF  (TTA»Z)  PAIRS 

no  o N e 1 »18,1 
READI2.inO)  TTA(N)»Z<N) 

100  FORMAT  (SX,F4,1»SX»F5.3) 

9 TTA(N)  a TTA(N) *3.141593/180. 


CALCULATING  R FOR  (R.TTA.Z)  COORDINATES 

no  10  N a 1.9.1 
10  R(N)  « 13.568  - .2516*Z(N) 

00  20  N ■ 10.18.1 
20  R(N)  ■ 16,625  - .4440*Z(N) 
no  30  L « 1.33.1 


RFADING  POINT  NUMBERS  FOR  END  POINTS  OF  LINE  SEGMENTS 

RFAn{2.101)  N1.N2 
101  FORMAT(5X.I3.5X.I3) 


CALCULATING  LENGTH  OF  LINE  SEGMENT 

30  ALONG(L)  » SORT ( ( R ( N2 ) *COS ( TTA ( N2 ) ) - R ( N 1 ) *COS ( TTA ( N1 ) ) 

1 ♦*?,  + ( R(N2)*SIN(TTA<N2)  ) - R<N1  )*SIN<TTA(,\1 ) ) ) **2. 

2 + ( Z(N2)  - Z(N1)  )**2.  ) 


DEFINING  (X.Y)  COORDINATES  OF  FIRST  TWO  POINTS 

X( 1 ) « 0. 

Y(1 ) • 3.057 
X<10)  - 0. 

Y(10)  - 0, 


CALCULATING  VERTEX  ANGLE  BETWEEN  SEGMENTS  . AND 

BESI  AVAIUBLE  COPY 
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C 

c 


c 

c 

c 

c 


(X,v)  CnORniMATES  OF  NFXT  POINT 


ALPHA  = 0,  I 

no  /.O  J ■ : 

r = ( (ALONG«?*J>  )**2,  “ (ALONG!?  *J+1))*«2.  - ( ALONG ( 2* J-1 ) ) ! 

1 **2.  ) / ( -2.*( AL0NG(2»J+1 ) )*AL0NG(2*J“1 ) ) I 

IF  (D 

•*4  V = + ATAN  ((SORTd.  - C»r))/C)  | 

GO  TO  “Xf  ■ 

■Xt  V = ATAN  ((SO<RT(].  - r*C))/C) 

•*A  i 

IF(M)  37.37,3fl  i 

N = ( J-^71  ) / ? ] 

ytrvH  s X(N-IO)  + {AL0NG(2*J+1)  )*SIN(V-ALPHA) 
v(N)  ® Y(N-IO)  - (AL0NG(2«J+1 ) ) »COS ( V-ALPHA ) 

GO  TO  40  • I 

•XA  N * (J+?)/? 

X(\)  a X(N+9)  + ( ALONG( 2*J+1 ) )*SIN(V-ALPHA)  ] 

Y(N)  a Y(N  + 9)  t-  ( ALONG!  ?*JTl  ))  «CnS  (V-ALPHA  ) j 

40  ALPHA  a V - ALPHA  I 


calculating  (X,Y)  COORDINATES  OF  POINTS  19  AND  20 

41  AR  a ATAN!  ! X ( Ifl )-X ! 9 ) ) / ! Y! 9 )“Y ( Ifl  ) ) ) 

X!70)  a X!19)  + ?,18R*SIN! 1.1529-AP) 

Y!?0)  a Y!1A)  + ?,18B*C0S( 1.1529-AP) 

X!1Q)  a X!20)  - 1.7114SIN!AP) 

V!lo)  s Y(?0)  + 1 ,711»C0S( AP ) 


dattfrn  ri^fifctfo  and  translated 

no  AG  N a 1,20»1 
4G  Y(N)  a A,  - Y(N) 


Plotting  pfgins 


CALI  SCALE  !1.0»1.0.0,»0.  ) 
CALI  FPLOT  !-?.X! 10) ♦Y! 10)  ) 
no  GO  N a 1,9*1 
CALI.  FPLOT!C»X!N)  ,Y(N)  ) 

AO  CALL.  POINT!  0) 

G1  CALL  FPLOT  ! 0 . X ! 1 9 ) ♦ Y ! 1 9 ) ) 
CALL  POINTIO) 

CALL  FPLOT  !0.X!20) #Y!20) ) 
CALL  POINTIO) 
no  70  I a 1,9,1 


BEST  AVAIUBLE  COPY 
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N'  * 19-1 

CALL  FPLOT(0»X(N) »Y(N)  ) 

70  CALL  POINT(O) 

CALL  PPL0T(0,X(1)  ,Y(1)  ) 
no  BO  I a l.lft*! 

0 a (-] . )**I 

IF(D)  71.71t7? 

71  N a 11  + ( I-l  ) /2 
00  TO  RO 

7?  N a (I  + ?)  /2 
BO  CALL  FPL0T(0»X(N)  ♦Y(N)  ) 
XX  * X(?0)  +2. 

CAll  FPLOT  (1»XX,0.) 

CALL  FXIT 

FNO 


BBt  AVAILABLE  COPY 


DATA 


N 

TTA 

(Degrees) 

Z 

(Inches) 

1 

0. 

0. 

? 

?. 

.17ft 

3 

4. 

.3P0 

4 

• 

.6?ft 

*) 

R. 

.920 

ft 
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